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AN INTERPOLATION FORMULA FOR "EQUIDIS- 
TANT" FREQUENCY DISTRIBUTIONS. 

By Harry Langman, Ocean Accident & Guarantee Corporation, 
New York City. 



A frequency distribution is conveniently represented by 
the area under some continuous curve. Moreover, the total 
frequency is usually included between two distinct values of 
the independent variable, outside of which there is no fre- 
quency. Hence, if x denote the independent variable, and 
V =/(*) denote the curve bounding the distribution, we must 
have 

(1) /(«) =/(&)= 0, 

where a and b are the extreme values of x. 

The frequency between any two values of x, a and 0, is 



given 



f 

by / y dx. In the usual case the values of this integral 



are known for each of the sub-intervals obtained by dividing 
the entire internal [a, b] into equal parts. Let m be the 
number of sub-intervals, and for convenience choose each 
sub-interval of unit length. Then 

(2) b = a+m. 

Let at, where fc = 1, 2, 3, •••, m, be the known frequencies over 
each of the sub-intervals [a, o+l], [a+1, a+2], ••; [a+k — 1, 
a+k], •••, [a+m— 1, a+m]. Then 

a + k 

(3) a*= ydx;k = l,2, 3, •■■, m. 

o + ik- 1 

The problem then is to interpolate frequencies for other 
than the equidistant values of the independent variable. 
This consists in finding some function f(x) which will satisfy 
the ra+2 conditions of (1) and (3).* Since this is infinitely 

*Cf. the problem* fiist solved by Lagrange, Oeuvres 1, p. 87. Also articles in Enzyk. d. Math. 
Wissen. IIA, 9a, p. 647. 
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indeterminate, the problem reduces to finding a form for 
f(x) which can be handled practically, and which can be applied 
in all cases. 

For this purpose we may choose for y = f(x) the function 

(4) y= 2 — -an'sm- - — ■ 

By (2), this evidently satisfies the conditions (1). The condi- 
tions (3) then become 

(5) ak= 2 Or- cos- — — cos ; « = 1, 2, 3, ••-, m. 

n — i L. to to _J 

The m conditions (5) are just sufficient to determine the m 
constants a n . The problem now is to evaluate these coef- 
ficients. 

Suppose the m equations of (5) to be written in succession 
and that we add the first k of them together, for all values 
of k. Then, if we put 

(6) s* = ai+02+a 3 + •• - +oi, 

n=m r ™_ I n=>m n=m ^_ 

sj= 2 a„- 1 — cos — k = So,- 2a a -cos — k. 

»-l L »» J n=l n-1 »» 

Letting 

(7) A = 2 a,, 

n=l 

the last may be written 

n=rti nT 

(8) Sk = A— 2 a n cos — ■ Jfc; fc = l,2, 3, •••, m. 

n=l m 

The to conditions of (8) replace those of (5). 
Let r be any integer from 1 to to. Then 

(9) r 2L m; n^m. 
Now consider the to quantities 

(10) 2 cos -k:k = l, 2, 3, •••, to. 

TO 

Multiplying each of the equations of (8) by the corresponding 
quantity of (10) and adding the resulting equations we obtain 

h—m y— n=m 

(11) 2 2 st cos- k<~A-I(r)+ 2a„-K(n,r), 

fc=l m n-l 
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where 



fc=M j. 



(12) f(f)=2 Scos-fc, 



h=\ 



m 



k ~ m n-ir , rir. 



and 

(13) JC(n,r)--2 2 cos — fc-cos-fc. 

h=\ m m 

From trigonometry, we have 

1 Sin(g0+|) 

(14) 2 cosp0=-±+ ^ fZ - 

" =1 J 2 sin | 

Putting q = m, and = — in (14), (12) becomes 

frh 



sin »+„- 
/(r , = _ 1+ _l 2mJ. 



(15) J(r) = -1 + 

sin^— 
2m 

By (9), the denominator of (15) never vanishes. Hence: 

(16) if r is even, I(r) =0; if r is odd, 1(f) = —2. 
From trigonometry, 

2 cos a;- cos y = cos (x+y)-\-cos (x—y). 
Applying this, (13) becomes 

(17) 2£ (w, r) = — 2 cos - — !-^-fc+cos — k 

= — 2 cos — k— 2 cos 1 — — k 

Again using (14), (17) becomes 



(18) K(n,r) = i- 



. f, , . . (n+r)x"| 
1 ml [( n+r ^ + -2^-J 



_ . {n+r)ir 
2 sin — - — — 
2m 

+1- 


sin 


(n r),r+ 2m 




2sin (w - ? - )7r 



m 
Now suppose n+r is even; then n— »• is even, and the angles 



77,— T 7l~\~T 

or -J; — is integral. By (9), — — is integral only when n=r 

71 — T 

= m; -„— is integral only when n=r. Hence we have: 

(19) if n+r is even, and n-\-r, K{nir) =0. 

Now suppose n+r is even and n=r. In this case (18) is 
indeterminate. However, on returning to (17) we obtain 



.. sin 2rir+ — 

1 I m I 1 

I = != =! — «i.= — 



. rir 
sin — 

(18a) if(r,r) = ^ *= l^J_ m= t TO . 

v v ' 2 . . rir 2 n . rir 

2 sin - 2 sin — 

m m 

Hence K= —m, unless the denominator is zero. This occurs 
only when r=m. Hence we have: 

(20) if r =1= m, K(r, r) = -m. 

Suppose n=r = m. Substituting in (17) we get at once: 

{21) K(m,m) = -2m. 

Consider n+r odd. Then n—r is odd, and the angles 
(n+r)ic and (n—r)ic in the brackets of (18) may each be 
replaced by t. It is clear, furthermore, that in this case the 
denominators of (18) cannot vanish. Hence we have 

(22) if n+r is odd, K(n, r) =2. 

For definiteness, we shall now suppose w odd. Then we 
may conveniently represent the results of (16), (19), (20), (21), 
(22); in the following table: 

n = r K=—m 



r even 



n even t 

' n=|=r K= 

n odd K= 2 



[1 = 



125] 



r odd 
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( r = m K=—2m 

n = r \ 
n odd \ [r=|=ra K=— m 

K= 
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\l=-2. 



n=\~r 



weven 



K = 



We note that for each value of r from 1 to m we get a set of 
multipliers in (10), and a new equation (11) from equations (8). 
We thus obtain m new equations to replace those of (8). 
That the new set is equivalent to (8) is shown by the fact that 
a definite solution is obtained. 

We may put 
(23) M = a 1 +a 3 +a !> +~-+a m : 

N =a 2 +a 4 +a 6 + •••+<*„,_!. 

Then 

A = M+N. 

Consider first the irregular case r = m. From the table we 
obtain: 1= — 2; for n even, K = 2; for n odd and n=|=m, K = 0; 
for n = r = m, K= —2m. Hence equation (11) takes the form 

2 S s*cosxfc= — 2A — 2mam+2(a i +a i +a 6 + ••■ +a m - 1 ) 



jb=i 



= -2A+2N -2ma m , 



or 



(24) 



s (-l)*«t= 



— M—mom- 



Let r assume any of the even values from 2 to m— 1. From 
the table, 7 = 0; K = 2 if n is odd; K = if n is even and n=\-r; 
K=—m if n=r. Hence (11) becomes 



h ** m rw 

2 2 St - cos — k- 

4=1 w 



■ 2(a!+a 3 + • • • +a m ) — ma, 



m- 



(25) =2M-ma r ;r = 2, 4, 6, 

Let r assume any of the odd values allowed by (9), exclusive 
of r=m, already considered. From the table we then obtain, 
7= —2; K = 2 if n is even; 7£ = if n is odd and n-\--r; K= —m 
if n=r. Hence (11) becomes 
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(26) 22 srcos — k = — 2A+2(«2+a 4 + •••+a m -i)— ma, 

= -2M-ma r ; r=l, 3, 5, •••, m-2. 
We have then m equations in (24), (25), (26). We now pro- 

ceed to determine M. Adding together the — — equations 

of (26), we obtain 

, m—l 

2i=m (Of 1 \ h 

2 2 2 .-ifcos ^— ; =-(f»-l)Jf-f«(oi+o, 



(=1 k=l 



m 



= — (m—l)M—m(M—a m ) 



(27) =-(2m-l)Jlf+OTa m 

Adding (24) and (27), we obtain 

1 



, »— l 

t r- 



^ M=- 2m 



k~m 2 k=m (Of i\ 1 

2(-l)*«+2 2 2 & • cos *- ; k . 

t-i (=1 *=i OT j 



This gives us M in terms of definitely known quantities. 
From (24), (25), (26), the coefficients a n are then determined. 
We can, however, decidedly simplify the expression for M in 
(28). 

From trigonometry, 

""' fn lNfl sin 20 

2 cos (2m- 1)9= - . ' 

tt=1 2 sin0 

unless sin0 = O. Using this, (28) becomes in order 



2m 



m—l 
k=m k=m 2 j 

2 (-l)*s*+2 2 st 2 cos (2<- 1)1* 

t=i *=i 1=1 m 

m—l 



t= 
1 fi—m k=m— 1 2 J. 

= ~o L 2(-l)*s*+2 2 « 2 cos(2/-l) — 

2m \k=\ fc=l 1=1 m 

m—l 



2 

+2 5 m 2 cos(2<-l)x 
(=1 
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. (m— 1)t4 
i (i-in h=m-l sm 

-=M 2 (-1)^+2 2 s k ±--( m -l) Sm 

2 « W i-1 2 on I* 



1 (i=ro 



7rfc , . irk 

k-sjn—l Bin7T«*cos cos 7r«*ein — 



fc=l 



k=l 



. irk 
sin — 









— (m 


- l)s m 


1 

1m 


2(-l)*5»- 2 (_i)*s t -(m-l>»[ 


1 
2m 


— 5m— ?W5 m +5mj 
J 


(29) = £.* 


From (24), (25), (26), we now obtain the coefficients a n in 


the following forms : 


1 f ? k=m 1 

«.= "£+ S(-D*s* ; 
™[2 j_, j 




, N 2.U m *-- rir, 
< 3 °) ar= ml2-\!/* C0S m fc 


, r=2, 4, 6, — , w-1; 




Or = 


2 


s m . *= m rir, 

7r+ 2 SjfcCOS — fc 


, r=l, 3, 5,-", m — 2. 





If in (4) a m be replaced by — , the three forms (30) for the 
coefficients a, can be expressed in the single form 



(31) 



r Sm £ 



ft«m 



rr 1 



Or=( — l) r 2« cos — k. 

m m k=1 m 



Suppose now m is even. Then we obtain in manner parallel 
to the case m odd : 



(28') 
(29') 



,. 1 J k Z m (2f — 1)tt t 

M= 2 2 Sk cos- —k 

m,_, ,, m 



Sm 

"2" 



*In practical computation, the two values found for M serve as an invaluable check. 
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The coefficients a n then become 
a« = - £- 2(-l)*5 t ; 

(30') o, = - £- 2 it cos-* ,r = 2,4, 6, ■■■,m-2; 
m [% 1=1 *» J 

a r = 7T+ 2 mcos — * , r=l, 3, 5, •••, w— 1. 

m\2 k=1 m J 

From (30') it is clear that (31) gives the coefficients a n also 
when m is even. Having expressed the coefficients a n in 
terms of the definitely known quantities s&, the function (4) is 
then completely determined.* 

In practical computation, the labor involved consists pri- 

marily in computing the values of 2 Sh cos — * for all values of 

i-i m 

tit 
rfLm. This requires the m 2 values of s* cos — *. But it is 

m 

fir 
readily observed that all the m 2 values of cos — * are given by 

m 

only of them (if m is odd). Furthermore, Sk is a "con- 

stant multiplier" for each of the values of cos — *, r^m. In 

m 

practice, further use can be made of the simple relations be- 
tween the trigonometric functions to materially shorten the 
work. The usual case arising for interpolation, is to obtain 
the frequencies for midway divisions of the sub-intervals. 
The computation made necessary for this is about equivalent 
to that in obtaining the coefficients a n from the given fre- 
quencies a*. 

*If the number of sub-intervals of [a, b], m, be indefinitely increased, it may be shown without 
difficulty, as is to be expected, that the function (4) becomes a Fourier series. 



